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The many-multiplet method applied to high redshift quasar absorption spectra has indicated a 
possible time variation of the fine structure constant. Alternatively, a constant value of a is consis- 
tent with the observational analysis if a non-solar isotopic ratio of 24,25,26 Mg occurs at high redshift. 
In particular, a higher abundance of the heavier isotopes 25,26 Mg are required to explain the ob- 
served multiplet splitting. We show that the synthesis of 25 ' 26 Mg at the base of the convective 
envelope in low-metallicity asymptotic giant branch stars, combined with a simple model of galactic 
chemical evolution, can produce the required isotopic ratios and is supported by recent observa- 
tions of high abundances of the neutron-rich Mg isotopes in metal-poor globular-cluster stars. We 
conclude that the present data based on high redshift quasar absorption spectra may be providing 
interesting information on the nucleosynthetic history of such systems, rather than a time variation 
of fundamental constants. 
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Over the last several years, there has been considerable 
excitement over the prospect that a time variation in the 
fine structure constant may have been observed in quasar 
absorption systems Q] - [H- While many observations 
have led to interesting limits on the temporal variation 
of a (see Iff for a recent review) , only the many-multiplet 
method [3 has led to a positive result, namely that ^ = 
(0.54 ± 0.12) x 10~ 5 over a redshift range of 0.5 < z < 3, 
where 5a is defined as the present value minus the past 
one. 

The sources of systematic errors in this method have 
been well documented 0,3 (see also Q)- Here, we would 
like to focus on one of these sources of systematic error 
for which there is recent evidence of a new interpreta- 
tion, namely the isotopic abundances of Mg assumed in 
the analysis. The analyses in - p|, have assumed ter- 
restrial ratios for the three Mg isotopes. They have also 
shown that had they neglected the presence of the neu- 
tron rich Mg isotopes, the case for a varying a would 
only be strengthened. They further argued, based upon 
the galactic chemical evolution studies available at that 
time, that the ratio of 25 ' 26 Mg/ 24 Mg is expected to de- 
crease at low metallicity making their result a robust and 
conservative one. 

In this paper, we will show that it is in fact quite plau- 
sible that the 25 ' 26 Mg/ 24 Mg ratio was sufficiently higher 
at low metallicity to account for the apparent variation 
in a. As such, we would argue that while the many- 
multiplet method of analysis does not unambiguously in- 
dicate a time variation in the fine structure constant, it 
can lead to important new insights with regard to the 
nucleosynthetic history of quasar absorption systems. 

Regarding the observations of Mg isotopes, Gay and 
Lambert |8| determined the Mg isotopic ratios in 20 stars 
in the metallicity range —1.8 < [Fe/H] < 0.0 with the aim 
of testing theoretical predictions [!j ■ (The notation [A/B] 



refers to the the log of the abundance ratio A/B rela- 
tive to the solar ratio.) Their results confirmed that the 
25 > 26 Mg abundances relative to 24 Mg appear to decrease 
at low metallicity for normal stars. Although many stars 
were found to have abundance ratios somewhat higher 
than predicted, even the 'peculiar' stars which show en- 
richments in 25 ' 26 Mg do not have abundance ratios sub- 
stantially above solar. 

Recently, however, a new study of Mg isotopic abun- 
dance in stars in the globular cluster NGC 6752 has been 
performed [lOj . This study looked at 20 bright red gi- 
ants which are all at a relatively low metallicity adopted 
at [Fe/H] = -1.62. These observations show a consider- 
able spread in the Mg isotopic ratios which range from 
24 Mg: 25 Mg: 26 Mg = 84:8:8 (slightly poor in the heavies) 
to 53:9:39 (greatly enriched in 26 Mgh The terrestrial 
value is 24 Mg: 25 Mg: 26 Mg = 79:10:11 (n|. Of the 20 stars 
observed, 15 of them show 24 Mg fractions of 78% or less 
(that is below solar), and 7 of them show fractions of 
70% or less with 4 of them in the range 53 - 67 %. This 
latter range is low enough to have a substantial effect 
on a determination of a in quasar absorption systems 
if the same ratios were to be found there. A previous 
study [13 also found unusually high abundances of the 
heavy Mg isotopes in M13 globular-cluster giants. Ra- 
tios of 24 Mg: 25 ' 26 Mg were found as low as 50:50 and even 
44:56. Similar results were very recently found in |13| . 
According to Q, raising the heavy isotope concentra- 
tion to 24 Mg: 25 ^ 26 Mg = 63:37 would sufficiently shift the 
multiplet wavelengths to eliminate the need for a vary- 
ing fine structure constant. While dispersion in the data 
could be a symptom of systematic errors often occurring 
when data from several samples are combined, real dis- 
persion is a signal that the observed abundances have 
been affected by local events. 

Available calculations of Type II supernova yields in- 
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put into chemical evolution models |9j, as well as obser- 
vations of the Mg isotopic abundances in relatively low 
metallicity stars, support the idea that the heavy Mg 
isotopes were rarer in the past. Nevertheless, this con- 
clusion is very sensitive to the the star formation history 
in object under consideration. Intermediate mass stars in 
their giant phase, are expected to be efficient producers 
of 25 ' 26 Mg [TiL Ha . fig and the recent data from a low- 
metallicity globular cluster 0, indicates substantial 
amounts of variation in the 25 ' 26 Mg/ 24 Mg ratios as well 
as a number stars highly polluted in the neutron rich Mg 
isotopes. These observations show isotopic ratios consid- 
erably higher than the ratios predicted in zero metallic- 
ity supernovae and they conclude that asymptotic giant 
branch (AGB) stars may be responsible for this contam- 
ination. 

Mg is produced in both Type I and Type II super- 
novae. In Type II supernovae, it is produced Il7j in the 
carbon and neon burning shells with an abundance some- 
what less than 10% of the oxygen abundance produced 
in massive stars. However, not much 25 ' 26 Mg is produced 
in conventional stellar evolution at low metallicity. This 
is because the isotopes 25 > 26 Mg are produced primarily 
in the outer carbon layer by the reactions 22 Ne(a,n) 25 Mg 
and 25 Mg(n,7) 26 Mg. According to the models of Woosley 
and Weaver , solar metallicity models produce final 
Mg isotope ratios reasonably close to solar values, while 
more massive stars tend to be slightly enhanced in the 
heavy isotopes, ( e.g., the 25 M Q model gives a ratio of 
65:15:20). Furthermore, the abundance of 25 > 26 Mg scales 
linearly with metallicity in the carbon shell. Hence, it 
would naively be expected that the ejecta from the first 
generation of supernovae would show a severe paucity of 
25 - 26 Mg. 

Much of the solar abundance of Mg is produced in 
Type la supernovae but with Mg to Fe ratios below solar. 
For ex amp le, the models of Thielemann, Nomoto and 
Yokoi 18], give [Mg/Fe] ~ -1.2. Due to the absence 
of free neutrons, essentially no 25 > 26 Mg are produced in 
Type la supernovae. 

The results of chemical evolution models tracing the 
Mg isotopes were presented in [fj. These models clearly 
show the effect of the low yields of 25 > 26 Mg at low metal- 
licity and predict [ 25 ' 26 Mg/ 24 Mg] < -0.8 at [Fe/H] < -1. 
The result of is essentially reproduced in the dashed 
curve shown in Figure 1 and discussed below. 

Based on the conventional theory of Mg production 
(e.g. [T3I, models of galactic chemical evolution j^] and 
previously available data on the isotopic abundance of 
Mg , the adoption of solar isotopic Mg ratios by Q - 
£| in the many-multiplet analysis would appear to be a 
safe and conservative one. These models indicate that 
the 24 Mg: 25 < 26 Mg ratio was higher than 79:21 in the past 
and it was shown Q that a higher ratio only strengthens 
the case for a varying a. These models, however, do not 
include contributions from intermediate mass stars as we 
now describe. Recall that increasing the abundances of 
the heavier Mg isotopes would yield a larger value for 



a, and a ratio of 24 Mg: 25 < 26 Mg = 63:37 is sufficient to 
obviate the need for a varying fine structure constant. 

Recently, it has been appreciated 0, 0, 0] that in- 
termediate mass stars of low metallicity can also be ef- 
ficient producers of the heavy Mg isotopes during the 
thermal-pulsing AGB phase. Heavy magnesium isotopes 
(and to some extent silicon isotopes as well) are synthe- 
sized via two mechanisms both of which are particularly 
robust in 2.5-6 Mq stars with low metallicity. Such low- 
metallicity objects, indeed are precisely the kinds of ob- 
jects which ought to produce the abundances observed in 
QSO narrow-line absorption systems at high redshift. 

One process is that of hot-bottom burning. This pro- 
cess is also believed to be a copious source of lithium 
in low metallicity stars (cf. 0|) During the AGB phase, 
stars develop an extended outer convective envelope. Ma- 
terial in this convective envelope is mixed downward to 
regions of high temperature at the base. Of particular 
interest for this paper is that the base of the envelope is 
more compact and of higher temperature in low metal- 
licity stars than in stars of solar composition. This can 
be traced to the decreased opacity of these objects. Fur- 
thermore, these stars would also have a shorter lifetime 
because they are hotter. Low to intermediate mass stars 
would contribute to the enrichment of the interstellar 
medium considerably sooner than their higher metallicity 
conterparts. 

Because these stars become sufficiently hot (T > 7 x 
10 7 K), proton capture processes in the Mg-Al cycle be- 
come effective. Proton capture on 24 Mg then leads to the 
production of 25 Mg (from the decay of 25 Al) and to 26 Al 
(which decays to 26 Mg). 

A second contributing process occurs deeper in the 
star during thermal pulses of the helium-burning shell. 
The helium shell experiences periodic thermonuclear run- 
aways when the ignition of the triple-alpha reaction oc- 
curs under electron-degenerate conditions. Due to elec- 
tron degeneracy, the star is unable to expand and cool. 
Hence, the temperature rapidly rises until the onset 
of convection to transport the energy away. During 
these thermal pulses, 22 Ne is produced by a captures 
on 14 N which itself is left over from the CNO cycle. 
Heavy magnesium isotopes are then produced via the 
22 Ne(ct,n)) 25 Mg and 22 Ne(a,7) 26 Mg reactions. It was ar- 
gued recently [15| , that in intermediate mass stars which 
experience a 3rd dredge-up, significantly greater amounts 
of 25 > 26 Mg are produced. A key point is that even though 
seed material is less plentiful in low metallicity stars, 
the reactions are very temperature sensitive. Hence, the 
increased temperature in the interior of low-metallicity 
stars more than compensates for the depleted seed ma- 
terial, leading to significant production of the heavy Mg 
isotopes. It has even been argued that these processes 
may also be net destroyers of 24 Mg 0, 0] due to the 
extreme temperatures attained. 

To illustrate the effects of producing enhanced abun- 
dances of 25 ' 26 Mg in intermediate mass stars, we show 
the results of a simple model of galactic chemical evo- 
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lution which traces the Mg isotopic abundances with 
and without the AGB source. When combined with 
the recent data showing such enhancements, our results 
suggest a plausible alternative for the interpretation of 
the quasar absorption system data based on the many- 
multiplct method. 

For our purposes a simple recalculation of the results of 
Timmes et al. [t| with and without the contribution from 
intermediate-mass AGB stars is sufficient. This allows us 
to make a direct comparison with the conclusions of the 
previous authors. 

The galactic chemical evolution model of Timmes et 
al. [2( is based upon exponential infall and a Schmidt 
star formation rate. We utilize a slightly modified model 
with updated yields [2(j which nevertheless reproduces 
the results of [2j in the appropriate limit. Hence, we 
write the evolution of the surface density a% of an isotope 
i as, 

da f 40 

— i = / B(t-T{m))V(m)X?(t-T(m))dm 
dt Jo. 8 

/•9.0 

+ / B(t-T{m))y{m)Xf GB (t-T(m))dm 

J2.5 

+mcoXl a Ri a - B{t)— + &i, gas , (1) 

a gas 

where B(t) is the stellar birthrate at t, Vf^m) is the initial 
mass function (IMF), Xf is the mass fraction of isotope 
i ejected from single star evolution and Type II super- 
novae, and r(m) is the lifetime of a star of mass m. The 
2nd term in Eq. ^ is the new contribution from AGB 
stars. For this purpose we adopt the AGB yields of [iH ]. 
The 3rd term in Eq. ^ is the contribution from type la 
supernovae where mco is the mass the exploding carbon- 
oxygen white dwarf, and Rj a is the SNIa supernova rate 
taken from ■ Finally, the last two terms represent the 
trapping of elements in new stars, and the galactic infall 
rate (presumed to be of primordial material) . 

As noted by a number of authors, the yields of heavy 
magnesium isotopes in AGB stars is extremely temper- 
ature sensitive, and hence rather sensitive to detailed 
physics of the stellar models. Moreover, there are reasons 
to expect that the initial mass function at low metallic- 
ity could be biased toward intermediate-mass stars. One 
argument for this is simply that with fewer metals, the 
cooling is less efficient in the protostellar cloud, so that 
a more massive cloud is required to form a star. To ac- 
count for these possibilities, we introduce a modest en- 
hancement of the IMF for intermediate mass stars. Such 
an enhancement has often been proposed and is moti- 
vated by models for star formation at low metallicity. For 
example, it has been invoked (2^ to account for a large 
population of white dwarfs as microlensing objects in the 
Galactic Halo. In fact, such a population of intermedi- 
ate mass stars was recently proposed .23] to explain the 
dispersion of D /H observed in quasar absorption systems 



[24| . Hence for the creation function, B^, we write: 

B(t)*(m) = Bi(t)*i(m) + B a (t)*a(m) (2) 
= B im - 235 + (5 2 /m)cxp(-log(m/5.0) 2 )/(2a 2 )). 

The IMF in Eq. J2J) accounts for a standard Salpeter 
distribution of stellar masses, 'J'i(m), with the addition 
of a lognormal component of stars peaked at 5 M Q , ^>2- 
The dimensionless width, a, is taken to be 0.07. This 
IMF is similar to a Gaussian with a width of 5a = 0.35 
Mq. For the normal stellar component we take the time 
dependence as 

Bt(t) = (1.0 - e-*/- 8Gyr ) - tot (t)[ - flO ./ -tatW] a , (3) 

while for the intermediate mass component we take 

B 2 (t) = b.he~ t l- 2G ^a tot {t)[a gas /a tot {t)] 2 . (4) 

This is very similar to the model used in |23|. It con- 
tains an early burst of intermediate mass stars peaked 
at 5 M which is exponentially suppressed after 0.2 Gyr. 
The second component describes standard quiescent star 
formation with a smooth transition from the burst. 

Figure 1 shows a comparison of our calculated magne- 
sium isotope ratio vs iron abundance. The solid curve 
shows the result of the model described above including 
the AGB contribution. The QSO absorption-line sys- 
tems in question have metallicities in the range from 
to -2.5 with a typical iron abundance of [Fe/H]~ —1.5. 
The mean isotopic ratio needed to account for the data of 
Q-H is 25 ' 26 Mg/ 24 Mg = 0.58 (shown by the solid hori- 
zontal line) with a la lower limit of 0.47 (dashed horizon- 
tal line). This figure clearly demonstrates that a plausi- 
ble model is possible in which a sufficient abundance of 
heavy Mg isotopes can be produced to both explain the 
observed globular-cluster data and the apparent trends in 
the many- mult iplet data or QSO absorption-line systems 
at high redshift. 

The behavior in the evolution of the heavy isotopes 
can be explained as follows: Initially, the production of 
25 > 26 Mg in the ejecta of intermediate mass stars is de- 
layed by their relatively long lifetime (compared to very 
massive stars). Initial contributions to the chemical en- 
richment of the interstellar medium comes from the most 
massive and shortest lived stars. In this model, the burst 
of intermediate mass stars begins to produce 25 ' 26 Mg at 
[Fe/H] > —2.5. At this stage, during the intermediate 
mass burst, 25 Mg and 26 Mg are copiously produced rela- 
tive to 24 Mg as per the yields of [Ty| ■ At higher metallic- 
ity, the ejecta from the standard population of (massive) 
stars which is poor in 25 > 26 Mg begins to dilute the ratio 
relative to 24 Mg, thereby producing the noticable bump 
in 25 - 26 Mg/ 24 Mg around [Fe/H] ~ -1.5. At late times, 
the effect of the early generation of intermediate mass 
stars is largely washed away. 

The dashed curve excludes the AGB yields and the 
intermediate mass component. It gives a result similar 
to that of We note that recently the AGB contribu- 
tion was included in a chemical evolution model |25| for a 
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FIG. 1: The chemical evolution of the 25 ' 26 Mg isotopes rela- 
tive to 24 Mg. The solid curve is our result based on eq. 
using the AGB Mg yields of [16J. The dashed curve is the 
result of turning off the AGB contribution and excluding the 
burst of intermediate mass stars. The horizontal lines indi- 
cate the ratio of 25,26 Mg/ 24 Mg necessary to explain the shifts 
seen in the many-multiplet analysis. 

normal stellar distribution. While the results showed sig- 
nificantly higher abundances of 25 - 26 Mg relative to 24 Mg 
than that given by the dashed curve, they were not high 
enough to account for the claimed variability in a. An 
enhanced early population of intermediate mass stars is 
therefore necessary. Nevertheless, this seems more plau- 



sible than a time varying fundamental constant. 

We have argued that previous models for the apparent 
broadening of the Mg multiplet in QSO absorption-line 
systems may have left out the important possible contri- 
bution from the production of heavy magnesium isotopes 
during the AGB phase of low-metallicity intermediate- 
mass stars. We have shown that a simple, plausible galac- 
tic chemical evolution model can be constructed which 
explains both the large abundances of heavy Mg isotopes 
observed in globular clusters and the large abundance 
necessary to explain the many-multiplet data. We note 
that the hot bottom burning process in AGB stars is 
also likely to have altered the Si isotopes as well by pro- 
ton captures on aluminum and silicon at the base of the 
convective envelope. Hence, it is possible that the sup- 
porting data from Si isotopes can be explained by this 
paradigm as well. Obviously more detailed work is war- 
ranted to clarify the ability of this mechanism to account 
for the data. Nevertheless, the model presented here is 
based upon plausible expectations of stellar and galactic 
evolution and should be taken seriously before demand- 
ing an alteration of any fundamental constant at high 
redshift. 
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